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Environmental estrogens have been shown to affect populations of aquatic 
organisms in devastating ways, including feminization of males, alterations in 
mating behaviors, and disruption of sexual selection. Studies have shown 17a- 
ethinylestradiol (EE2) exposure to induce female-like secondary sexual traits in 
male Gulf pipefish, changing how females perceive affected males. We aimed to 
understand the effects of EE2 exposure on the sex-role-reversed mating system 
and the strength of selection in Gulf pipefish. We used artificial Gulf pipefish 
breeding aggregations and microsatellite-based parentage analysis to determine 
maternity. We then calculated the opportunity for selection and selection differ- 
entials on body size for both sexes during three consecutive episodes of selection. 
Exposure to EE2 did not affect the strength of selection, likely due to the unusual 
sex-role-reversed mating system found in this species. With respect to multiply 
mated females, EE2 -exposed females produced more eggs with higher embryo 
survivorship than nonexposed females. Thus, short-term exposure to low con- 
centrations (2.0 ng/L) of EE2 in Gulf pipefish enhanced female reproductive suc- 
cess. However, higher EE2 concentrations (5.0 ng/L) caused complete 
reproductive failure in Gulf pipefish males. These results call for more work on 
the long-term effects of EE2 exposure in Gulf pipefish in artificial and natural 
populations. 
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Introduction 

Endocrine disrupting chemicals can mimic natural 
hormones and alter bodily processes regulated by the 
endocrine system, causing detrimental effects on reproduc- 
tion and hormone production (Orlando and Guillette 
2007). The earliest studies on the impacts of endocrine 
disruptors focused primarily on the negative physiological 
effects of these compounds on the reproductive systems of 
exposed animals, particularly in aquatic organisms such as 
frogs and fish (Allen et al. 1999; Iguchi et al. 2001). These 
studies documented partial feminization or complete sex 
reversal in exposed males and termination of egg produc- 
tion in females (Lange et al. 2001; Islinger et al. 2003; Xu 
et al. 2008). One endocrine disruptor that has received 
recent attention is a synthetic estrogen found in hormonal 
contraceptives called 17a-ethinylestradiol, also known as 



EE2. This contaminant is resistant to degradation in the 
body and ultimately ends up in the aquatic environment by 
passing through domestic wastewater treatment facilities 
and being released as a biologically active molecule in the 
effluent where it can accumulate to levels that cause delete- 
rious effects on exposed organisms (Kolpin et al. 2002). 
For instance, EE2 has been detected in US rivers at levels as 
high as 820 ng/L, European locations at 35 ng/L, and at 
42 ng/L in Canadian sewage treatment effluent (Ternes 
et al. 1999; Kolpin et al. 2002; Pojana et al. 2007). These 
high contamination levels raise questions regarding the 
types of effects that varying levels of EE2 exposure could 
have on populations of economically or ecologically impor- 
tant species that occupy such sites during all or part of their 
life cycles (Segner et al. 2003; Soares et al. 2009). 

Lower ranges of EE2 have been detected in surface waters 
near sewage treatment plants, generally ranging from below 
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0.5 to 5 ng/L in locations worldwide (Allen et al. 1999; 
Ternes et al. 1999; Johnson et al. 2000). Chronic exposure 
to lower concentrations of EE2 has been shown to cause 
entire populations of fathead minnow (Pimephales prom- 
elas) to stop reproducing after a single generation of 
5-6 ng/L exposure (Kidd et al. 2007). Relatively low levels 
of EE2, ranging from 0.5 to 5.0 ng/L, have been shown to 
have drastic effects on gene expression levels, reproductive 
development, and behavior in exposed populations (Nash 
et al. 2004; Larsen et al. 2008; Ferreira et al. 2009; Soares 
et al. 2009). Even if endocrine disruptors have no obvious 
morphological or physiological effects, changes in either 
gene expression levels or behavior may hinder the organ- 
ism's reproductive success. 

Exposure to endocrine disruptors, such as EE2, can alter 
an organism's mating behaviors at several different levels of 
biological organization. For example, endocrine disruptors 
can improperly induce or eliminate secondary sex traits, 
decrease sexual behaviors such as courtship displays or 
male territory aggression, disrupt communication between 
the sexes, or simply sterilize the exposed animals by inter- 
fering with the normal development of reproductive struc- 
tures (Orlando and Guillette 2007; Coe et al. 2008; 
Munakata and Kobayashi 2009; Saaristo et al. 2009b, 
2010). In medaka fish (Oryzias latipes), for instance, EE2 
has been shown to disrupt mating rituals, such as dancing 
and following with their mates; however, normal behaviors 
such as resting and swimming were not affected (Oshima 
et al. 2003). Endocrine disruptors have been shown to alter 
courtship displays and aggression in several other species of 
fish including goldfish (Carassius auratus), zebrafish (Danio 
rerio), three-spined sticklebacks (Gasterosteus aculeatus), 
salmon (Oncorhynchus masou), and guppies (Poecilia retic- 
ulata) (Bell 2001; Bjerselius et al. 2001; Kristensen et al. 
2005; Colman et al. 2009; Sarria et al. 2011a). Several of 
these studies demonstrated drastic changes in male behav- 
iors, such as males completely failing to court females or 
no longer defending their territories, at EE2 concentrations 
as low as 0.5-3 ng/L (Bjerselius et al. 2001; Oshima et al. 
2003; Colman et al. 2009). Sarristo et al. found that EE2 
exposure affects courtship, aggression, and parental care in 
male sand gobies, Pomatoschistus minutus (Saaristo et al. 
2009a, 2010). Additionally, male sand gobies exposed to 
EE2 experienced difficulty defending their nests, causing 
females to prefer nonexposed males when given the choice 
(Saaristo et al. 2009b). This study demonstrates the strong 
impact EE2 can have on an exposed male's mating success 
by disrupting the mating system and decreasing the impor- 
tance of secondary sexual characteristics in males. 

We chose to study the effects of EE2 on sexual selection 
in the sex-role-reversed Gulf pipefish, Syngnathus scovelli, 
because this species resides in areas that are likely to be 
affected by EE2 and has previously been the focus of many 
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behavioral studies. Gulf pipefish are found in seagrass beds 
along the Gulf of Mexico and Atlantic coastlines of North 
America (Dawson 1982) and depend on the seagrass com- 
munity for their habitat, food, and protection. This reliance 
on the seagrass ecosystem ties the organisms to locations 
near the coast, often in the vicinity of sewage treatment 
plants and other sources of environmental contamination. 
There are numerous benefits to using the Gulf pipefish 
mating system to test the effects of EE2 on the strength of 
sexual selection. The Gulf pipefish, like the other members 
of the family Syngnathidae, exhibits an evolutionarily novel 
trait in male pregnancy and as a result is sex-role reversed, 
meaning that females compete for access to mates and 
males act as the choosy sex during courtship. Females of 
this species are typically larger than males at sexual matu- 
rity, and males prefer larger females (Jones and Avise 
2001). Gulf pipefish are sexually dimorphic with respect to 
secondary sexual traits as well, with sexually mature females 
exhibiting iridescent bands on their trunk, a deeply keeled 
abdomen, and an enlarged, darkened dorsal fin (Dawson 
1982). None of these traits normally occur in males. The 
Gulf pipefish has a polyandrous mating system, where 
males typically mate with a single female per pregnancy but 
successful females transfer eggs to multiple mates (Jones 
and Avise 1997). Consequently, sexual selection acts more 
strongly on females (Jones et al. 2001). Pregnant male Gulf 
pipefish carry eggs in a sealed brood pouch, located on the 
ventral surface of his body, where the eggs that have been 
transferred from the female are encased and fertilized, 
resulting in assured paternity of the offspring (Jones et al. 
2001). Unsuccessful eggs, or those that do not develop, per- 
sist within the brood pouch until the end of the pregnancy, 
making it easy to accurately count the number of successful 
and unsuccessful eggs transferred (Paczolt and Jones 2010). 

Our work builds on several studies that have addressed 
the effects of endocrine disruptors on various species of 
pipefish (Ripley and Foran 2008; Sarria et al. 2011b). For 
example, Sarria et al. (2011c) found that the juvenile black 
striped pipefish, S. abaster, when exposed to both tributyl- 
tin (TBT) and EE2, altered their swimming patterns and 
their response to the mosquitofish (Gambusia affinis), a 
potential predator. The Gulf pipefish has recently emerged 
as a useful model for testing the effects of endocrine disrup- 
tors. For instance, EE2 levels of 1 and 100 ng/L have been 
shown to affect male morphology, gene expression levels, 
and mating dynamics in the Gulf pipefish (Ueda et al. 
2005; Partridge et al. 2010). As a result of EE2 exposure, 
male pipefish experience feminization by developing a 
deeper abdomen and iridescent bars, traits typically only 
seen in females. Female Gulf pipefish prefer nonexposed 
over exposed males, and exposed males were less likely to 
become impregnated when chosen as a mate (Partridge 
et al. 2010). These results from binary mate choice assays 
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indicate that mating dynamics in Gulf pipefish are affected 
by EE2 exposure. Sex-role-reversed pipefish have also been 
shown to have high levels of natural estrogen in brooding 
males, a reversal of the normal pattern of estrogen in male 
teleosts, suggesting exposure to a synthetic estrogen, EE2, 
might affect selection acting on male pregnancy (Mayer 
et al. 1993). The next question, which we address here, 
concerns how this disruption in mating preferences alters 
the mating system and the strength of selection in pipefish 
breeding aggregations. We set out to determine the effects 
of low levels of EE2 exposure on both pre- and postmating 
episodes of selection in artificial breeding colonies of Gulf 
pipefish in a laboratory setting. To accomplish this goal, we 
measured mating success, reproductive success, and 
embryo survivorship within small breeding aggregations of 
Gulf pipefish, which were either exposed or not exposed to 
low concentrations of EE2 and used these data to measure 
the effects of EE2 on the intensity of selection and other 
reproductive attributes of this sex-role-reversed species. 

Methods 

Gulf pipefish were collected from coastal seagrass beds near 
Aransas Pass, Texas (27°53'39.07"N, 97°7'51.69"W) from 
July through October 2010. Sexually mature males and 
females were separated by sex, acclimated to 26-ppt salinity 
tanks, and group housed in a flow-through system at Texas 
A&M University. All males were collected pregnant to con- 
firm a history of successful reproduction and were allowed 
to give birth in the laboratory. Males were used in the 
experiment within a month of their collection to ensure a 
recent pregnancy. The EE2 concentrations at the collecting 
site are currently unknown. However, the location was cho- 
sen because of its long distance from any outflows from 
sewage treatment plants. 

To determine the levels of EE2 exposure for our experi- 
ment, we performed a pilot study at EE2 concentrations of 
2 and 5 ng/L with a total of eight males and eight females 
per treatment along with a parallel control set of fish. We 
placed males into the experiment while still pregnant with 
their broods from the field and monitored their abilities to 
carry their broods to term and to become pregnant with 
subsequent broods. While control males gave birth and 
mated soon thereafter as anticipated, several males exposed 
to 5 ng/L of EE2 experienced difficulties in giving birth and 
had dead offspring in their pouches after the first few days 
of exposure. None of the males exposed to 5 ng/L of EE2 
had a second pregnancy in the laboratory and they 
appeared to resorb their brood pouches, apparently termi- 
nating their reproductive activities. Females, in contrast, 
continued to show courtship displays, including dancing 
and twitching, similar to the control females. Thus, we 
concluded from this pilot study that a concentration of 
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5 ng/L of EE2 would result in complete reproductive fail- 
ure for Gulf pipefish, apparently mediated by a loss of 
reproductive ability for males and chose 2 ng/L as our EE2 
concentration to allow for more typical pipefish mating 
and offspring development. 

Over the course of the experiment, we conducted seven 
experimental replicates and seven control replicates. Each 
replicate began with eight nonpregnant, adult males and 
eight adult females in a 100-liter tank. For the experimental 
tanks, we used a 2 ng/L EE2 concentration, whereas the 
control tanks were EE2 free. The 17a-ethinylestradiol pow- 
der, of 98% purity, was obtained from Sigma (St. Louis, 
MO, USA; # 028K1411, MW 296.4, CAS 57-63-6) and dis- 
solved in ethanol. Tanks treated with EE2 were initially 
dosed on the first day to obtain a concentration of 2 ng/L, 
and 10% water changes were conducted daily to maintain a 
constant level of 2 ng/L exposure as established by Par- 
tridge et al. (2010). Control tanks were treated with the 
same volume of ethanol (100 /(L) without 17a-ethinylestra- 
diol. We used 100-liter tanks that were optimized for 
length (92 x 30 cm), rather than height (38 cm), and were 
taller than previous tanks (28 cm) used for mating trials in 
Gulf pipefish (Paczolt and Jones 2010) to allow a large 
amount of surface area for fish matings while minimizing 
the amount of EE2 wastewater. 

On the first day of the experiment, fish were marked with 
three visible implant fluorescent elastomer tags (VIFE; 
Northwest Marine Technology, Inc., Shaw Island, WA, 
USA) for identification using the protocol of Wood and 
Martin-Smith (2004). After anesthetization with diluted 
clove oil, each fish was marked on its tail with at least one 
blue and one yellow band to minimize color differences 
among marks. The marking procedure produced no mor- 
talities, and we found no evidence for preferences for a sin- 
gle marking pattern during the experiment across both 
treatments (C: F 7M = 0.911, P = 0.5062; E: F 7>55 = 1.369, 
P = 0.2401). Fish were then randomly placed into treat- 
ments and were not significantly different in size across the 
two treatments (f-tests, males: n = 109, P = 0.6781; 
females: n=\U,p = 0.0904). 

For 3 days before the establishment of mixed-sex breed- 
ing aggregations in the 100-L tanks, we housed males and 
females separately from one another and exposed them to 
the desired level of EE2 (i.e., 2 ng/L for experimental 
animals and 0 ng/L for control animals). On the fourth day 
of the experiment, sexes were combined in 100-liter tanks 
and the EE2 treatments were maintained throughout the 
rest of the experiment. Each replicate with eight males and 
eight females per tank was given 2 weeks during which 
mating took place. Males were checked daily for pregnan- 
cies and sacrificed using MS222 on day eight of their preg- 
nancy. Eggs were dissected out of the male's pouch to 
determine the number of developing offspring and failed 
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eggs. The proportion of normally developing embryos in 
each brood is used as a measure of embryo survivorship. 
Four offspring were removed from the top and bottom of 
the pouch and preserved in ethanol to use for assigning 
maternity. All nonpregnant males and females were sacri- 
ficed on day 18 of the experiment. Dorsal fins from the 
adult fish were preserved in ethanol for DNA extractions 
using the Genomic DNA Purification Kit from Gentra sys- 
tems (Qiagen, Germantown, MD, USA). DNA was 
extracted from the embryos using a Chelex/proteinase K 
(20 mg/mL) extraction method. 

We conducted a microsatellite-based parentage analysis 
to determine maternity of the broods using three highly var- 
iable microsatellites (micro25.10, micro25.22, and micro22.3) 
previously developed by for S. scovelli (Jones and Avise 
1997). Paternity was already known because the offspring 
were dissected out of the male's brood pouch, where fertil- 
ization occurs and males are assured paternity (Jones and 
Avise 2001). Microsatellites were amplified for all adults and 
eight offspring per male with the exception of the eggs from 
males with zero offspring survivorship, and PCR products 
were sent for fragment analysis at the Cornell Life Sciences 
Core Laboratories Center. Microsatellite fragment sizes were 
measured by an Applied BioSystems 3730x1 DNA Analyzer 
and analyzed using Peak Scanner software (Life Technolo- 
gies, Grand Island, NY, USA). For each male's brood, a 
maximum of four alleles were represented, two of which 
matched the paternal genotype and two that represented the 
maternal genotype. Within each replicate, one of the eight 
females from the appropriate replicate was unambiguously 
matched with each set of embryos by exclusion. Maternity 
was easily assigned for all male broods, except for two males 
in the control tanks and one in the EE2 whose embryos 
failed to develop and thus were not amenable to microsatel- 
lite analysis. These three failed broods were excluded from 
further analysis. 

Fish were photographed on the first, eighth, and four- 
teenth days of the experiment to measure standard body 
length and depth using ImageJ (NIH, Bethesda, MD, USA). 
Standard body length was measured from the snout of the 
fish to the end of the caudal peduncle, and standard depth 
was measured from the anterior end of the dorsal fin to the 
base of the fish's ventral surface. We calculated the covari- 
ance between standard length and relative fitness to mea- 
sure the absolute selection differential, s (reported in cm), 
and standardized selection differential, s' (reported in units 
of phenotypic standard deviations; Lande and Arnold 
1983). Absolute and standardized selection differentials 
were calculated across three episodes of sexual selection, 
broken down into one premating episode pertaining to 
mating success and two postmating episodes, including eggs 
transferred per mate and embryo survivorship (Arnold and 
Wade 1984a). These latter two episodes could be considered 



either sexual or natural selection, depending on the mecha- 
nisms involved. We also decomposed the variance in rela- 
tive fitness, also known as the opportunity for selection (f), 
into components arising from these three episodes of selec- 
tion (Arnold and Wade 1984a,b). We calculated s, s', and I 
for each tank individually and reported the means and stan- 
dard errors for each episode across replicates. Relative fit- 
ness for each episode of selection was calculated separately 
for each replicate tank by dividing each measure of absolute 
fitness (i.e., number of mates, number of eggs transferred, 
or number of surviving offspring) by the corresponding 
mean absolute fitness across all the individuals of the same 
sex in the replicate. A more detailed analysis of the control 
tanks, including a more extensive discussion of the inter- 
pretation of the various phases of selection, has been pub- 
lished elsewhere (Rose et al. 2013). All other statistical tests 
were conducted using JMP 9.0 (SAS, Cary, NY, USA). 

Results 

Did EE2 have an effect on mating success? 
The presence of low levels of EE2 did not have an effect on 
the number of males that mated successfully. Each treat- 
ment saw a similar number of males become pregnant: 
91% of males mated in the control tanks with only five of 
the 56 males failing to become pregnant and 89% of the 
males mated in the EE2 treatment, where only six of 56 
males remained unmated (Fig. 1). Parentage analysis of 
eight eggs per male, four from each end of the brood 
pouch, confirmed that all males mated with a single female, 
regardless of treatment. Using an anova with tank as a 
random effect, we found no effect of treatment on male 
size. However, despite a small sample size of unmated 
males, we did see a significant difference in size of mated 
and unmated males across the entire experiment (treat- 
ment: -Fi,io9 = 0.005, P = 0.945; mating category: 
F 1,109 = 4.06, P = 0.047; treatment*mating category: 
F ljl0 9 = 0.074, P = 0.787). Mated males averaged 8.73 cm 
(SE: 0.15) in length, while unmated males were on average 
8.28 cm (SE: 0.25) in length. 

The number of mated females was slightly higher in the 
EE2 treatment than that in the control, with 61% of 
females mating successfully (34 of 56) in the EE2 tanks as 
compared to 53% (29 of 55) in the control (Fig. 1). The 
control had more multiply mated females (n = 14) than 
the EE2 replicates (n = 11), whereas the EE2 treatment had 
more singly mated females (n = 23) than the control 
replicates (n = 15). Using an anova with replicate as a ran- 
dom effect, we found no significant effect of treatment and 
whether or not the female mated on female body length 
(ancova: treatment: Fijio = 1.43, P = 0.254; female mating 
category: Fi >110 = 3.225, P = 0.074; treatment* female mat- 
ing category f ul0 = 0.448, P = 0.505). 
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Number of mates 

Figure 1 Histograms displaying mating success for males and females 
across both treatments with black bars for the control and gray bars for 
the EE2 treatment. The y-axis represents the frequencies, and the x-axis 
represents the number of mates. 

Did EE2 have an effect on reproductive success? 

To examine female reproductive success, we compared the 
number of eggs females transferred throughout the whole 
experiment as a function of treatment and whether the 
female was unmated, singly mated, or multiply mated. We 
also used female body length as a covariate, included 
replicate as a random effect, and retained only the 
significant interactions in the final model (ancova: treat- 
ment: J 7 !,!! = 2.17, P = 0.164; female mating category: 
^2,110 = 246.5, P < 0.0001; female body length: 
P ljll0 = 3.61, P = 0.06; treatment*females mating category: 
^2,110 = 4.48, P = 0.014). As expected, multiply mated 
females transferred a greater number of eggs to males over 
the course of the experiment simply because they had 
gained access to more total brood pouch space (female 
mating category: £2,110 = 246.5, P < 0.0001). On average 
across the treatments, multiply mated females transferred a 
total of 69 eggs compared with 32 eggs transferred per 
singly mated female. We also found a significant interac- 
tion between treatment and female mating category 
(treatment*female mating category: f 2,110 = 4.48, 
P = 0.014). This interaction was driven by the multiply 
mating females in the EE2 experiment transferring a 
significantly larger number of eggs than any other category 



of female, including multiply mating females in the control 
tanks (Fig. 2A). Thus, we did not see a difference between 
treatments in terms of the numbers of eggs transferred by 
singly mated females (Tukey's post hoc: n = 38, P = 0.91), 
but we did observe that multiply mated females in the EE2 
treatment transferred 15 more total eggs on average 
compared with multiply mated control females (Tukey's 
post hoc: n = 25, P = 0.05). 

We compared reproductive success on a per-brood basis 
across treatments and female mating categories, including 
singly or multiply mated. We used an ancova to examine 
the effects of treatment and female mating category on the 
number of surviving offspring with the total number of 
eggs initially transferred as a covariate. We included 
replicate as a random effect and reported only the signifi- 
cant interactions. In the control, females mating with 
multiple males experienced lower survivorship per embryo 
relative to singly mated females, whereas in the EE2 treat- 
ment, no such pattern was evident (ancova: treatment: 
P i,96 = 0.029, P = 0.87; female mating category: 
F h96 = 0.774, P = 0.38; total eggs transferred: F h96 = 1488, 
P < 0.0001; treatment*females mating category: 
F 1>96 = 7.66, P = 0.007). Thus, in control tanks, females 
appear to experience a trade-off between number of eggs 
transferred and embryo survivorship, but this trade-off 
disappears in the EE2 treatment, despite the fact that 
multiply mating EE2 females actually transferred a greater 
number of eggs than multiply mating control females 
(Fig. 2A,B). This pattern is also evident from an embryo 
survivorship standpoint (Fig. 2B): an anova shows a signif- 
icant difference in embryo survivorship across the 




Control EE2 Control EEZ 



Figure 2 Number of eggs transferred and proportion of surviving 
embryos for females in the EE2 and control treatments. The graph on 
the left (A) shows the number of eggs females transferred to their 
mates over the entire experiment. Solid bars represent singly mated 
females, and striped bars represent multiply mated females. The graph 
on the right (B) shows the proportion of surviving embryos for singly 
mated (solid bars) and multiply mated (striped bars) males in the control 
and EE2 treatments. Within each figure, bars with shared letters are not 
significantly different from one another (Tukey's post hoc test). Error 
bars represent one standard error from the mean. 
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experiment (-F 3> 9 5 = 3.59, P = 0.017). This pattern is a result 
of a significant difference between embryo survivorship for 
males mated with singly versus multiply mated females in 
control tanks (Tukey's post hoc: n = 47, P = 0.033), but 
not between males mated with singly and multiply mated 
females in EE2 tanks (Tukey's post hoc: n = 50, P = 0.446). 

Did EE2 have an effect on selection? 

The decomposition of selection differentials is presented in 
Table 1. For both treatments, the largest contribution to 
the total selection on body size in females was from the first 
episode of selection, mating success. Selection on body size 
resulting from mating success in the EE2 females was sig- 
nificantly positive at a. = 0.05, and we observed a similar, 
nonsignificant trend in control females. The selection dif- 
ferentials also provide evidence of a trade-off between pre- 
and postmating episodes of selection in control females, a 
result that was also apparent (and statistically significant) 
from our analysis of the relationship between multiple mat- 
ing and embryo survivorship. For control females, both 
episodes of postmating selection, resulting from eggs per 
mate and offspring survivorship, are negative and oppose 
the strong positive selection on body length from mating 
success (Table 1). Interestingly, we saw no evidence for this 
trade-off in the EE2 tanks, where selection differentials for 
premating sexual selection as well as for both postmating 
episodes of selection were all positive, favoring larger 



females (Table 1), but these results also must be interpreted 
with caution as the 95% confidence intervals for the selec- 
tion differentials for postmating episodes overlapped zero. 
Regardless, our results do clearly show that sexual selection 
on females was not dramatically reduced as a result of EE2 
exposure in Gulf pipefish. 

Patterns of selection in males differed considerably from 
those in females. For instance, no single episode of selec- 
tion made up the vast majority of the male selection differ- 
entials (Table 1). For both control and EE2 males, we see 
more of a balance between the contributions of the first 
two episodes. The second episode of selection measures 
reproductive success per pregnancy and suggests that larger 
males are receiving a greater number of eggs. This episode 
is statistically significant in control and EE2 treatments for 
the standardized selection (s') differential on male standard 
length and is probably best described as fecundity selection, 
because larger males, which typically have larger brood 
pouches, receive more eggs per mating (control: n = 49, 
r = 0.38, P = 0.007; EE2: n = 50, r = 0.29, P = 0.044). 

We also calculated the total opportunity for selection 
and found that the decomposed episodes of I paralleled the 
results from our selection differentials (Table 2). The total 
opportunity for selection, as well as I for each of the three 
individual episodes, did not differ across the treatments for 
either sex. However, similar to our total selection differen- 
tials, we see that in both the control and EE2 treatments, 
the total opportunity for selection is statistically greater in 



Table 1. Selection differentials broken down into three episodes of selection for the control and EE2 treatments. We calculated the absolute selec- 
tion differential (s), in cm, and the standardized selection differential (s'), in units of phenotypic standard deviations, for each tank separately and cal- 
culated means and 95% confidence intervals (bracketed values below the mean) across tanks. Episodes of selection include mating success, number 
of eggs transferred per mate, and embryo survivorship. 



Selection episode 



Males 
[95% CI] 



Male s' 
% [95% CI] 



Female s 
% [95% CI] 



Female s' 
% [95% CI] % 



Control 

Premating selection 

(mating success) 
Postmating selection 

(eggs per mate) 
Postmating selection 

(embryo survivorship) 
Total selection differential 

EE2 

Premating selection 

(mating success) 
Postmating selection 

(eggs per mate) 
Postmating selection 

(embryo survivorship) 
Total selection differential 



0.034 41.9 

[-0.012, 0.079] 

0.037 46.0 

[-0.001, 0.075] 

0.010 12.1 

[-0.044, 0.064] 

0.080 100 

[-0.002, 0.163] 

0.037 26.8 
[-0.006, 0.081] 
0.100 71.4 
[-0.016, 0.215] 
0.002 1.8 
[-0.015, 0.019] 
0.139 100 
[0.013, 0.266] 



0.087 

[-0.046, 0.220] 
0.051 

[0.006, 0.096] 
0.024 

[-0.053, 0.102] 
0.162 

[0.004, 0.321] 
0.082 

[-0.053, 0.216] 
0.133 

[0.005, 0.261] 
0.0 

[-0.019, 0.019] 
0.215 

[0.067, 0.363] 



53.5 0.140 

[-0.015, 0.295] 
31.4 -0.026 

[-0.093, 0.041] 
15.0 -0.016 

[-0.029, -0.004] 
100 0.097 

[-0.046, 0.241] 

38.0 0.120 

[0.029, 0.211] 
62.0 0.029 

[-0.022, 0.080] 
0.0 0.006 

[-0.009, 0.021] 
100 0.155 

[0.037, 0.272] 



143.5 0.209 155.1 

[-0.011, 0.429] 

-26.9 -0.044 -32.8 

[-0.174, 0.086] 

-16.6 -0.030 -22.4 

[-0.056, -0.005] 

100 0.135 100 

[-0.072, 0.341] 

77.4 0.236 86.1 

[0.037, 0.436] 

18.7 0.030 10.8 
[-0.050, 0.109] 

3.9 0.008 3.1 
[-0.014, 0.031] 

100 0.274 100 
[0.050, 0.499] 
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females than that in males as evidenced by nonoverlapping 
95% confidence intervals (Table 2). When we break down 
the sources of variation responsible for / into episodes, 
there are very similar patterns across each treatment within 
the sexes when compared to our decomposed selection dif- 
ferentials. Similar to s and s', we find that the first episode 
of selection, resulting from variance in mating success, is 
responsible for the majority of variance in female fitness 
and is statistically significant in both the control and EE2 
replicates. On the other hand, the second and third epi- 
sodes of selection only represent 4-4.6% and 0.6-0.9%, 
respectively, of the total opportunity for selection in 
females across treatments (Table 2). The total opportunity 
for selection in males was significantly lower than that in 
females, and the decomposition in males reveals different 
contributions of the various episodes of selection. In gen- 
eral, we found similar patterns in control and EE2 treat- 
ments for males. In particular, mating success (Ii) and 
number of eggs per mate (I 2 ) make the largest contribution 
to variation in fitness for males. Embryo survivorship (7 3 ) 
made a small (i.e., 4-17.5%) but statistically significant 
contribution in both treatments (Table 2). 

Our final metric related to sexual selection was the Bat- 
eman gradient, which we calculated for females and com- 
pared across treatments (Fig. 3). Our results show that the 
Bateman gradient for females is significantly steeper in the 
EE2 treatment than it is in the control treatment (ancova 
with replicate as a random effect: P = 0.009), implying that 
sexual selection may be slightly stronger in the EE2 treat- 
ment than in the control treatment. 




♦ Control 
□ EE2 

Linear (Control) 

Linear (EE2) 



-20 



1 2 3 

Mating success 



Figure 3 Absolute Bateman gradients for females in the EE2 and 
control treatments. The black line represents the linear regression of 
total number of offspring on number of mates for the control 
treatment, whereas the gray line shows the same linear regression for 
the EE2 treatment. The Bateman gradient for females in the EE2 
treatment is significantly steeper than the Bateman gradient for females 
in the control replicates (ancova, with replicate as a random effect: 
P = 0.009). 



Discussion 

Exposure to EE2 has previously been shown to alter the 
mating behaviors of Gulf pipefish in binary choice tests, 
which suggested that EE2 could potentially disrupt the 
mating system and alter sexual selection. While we did find 
several changes in the mating system of Gulf pipefish as a 
result of low levels of EE2 exposure, we did not see a signif- 
icant change in the strength of selection for either sex. The 
presence of 2 ng/L of EE2 in our experimental treatments 
did not hinder the ability of males to become pregnant. 
However, exposure to higher concentrations of EE2, 
including 5 ng/L in our pilot study and 100 ng/L in previ- 
ous studies, even for a short period of time, significantly 
impacted male reproductive potential (Partridge et al. 
2010). While EE2 concentrations of 5 ng/L are typically 
considered to be on the lower end of the range of EE2 levels 
detected in the natural environment, studies have shown 
that chronic exposure to these relatively lower levels of EE2 
can sterilize male fishes (Kolpin et al. 2002; Kidd et al. 
2007). Thus, high levels of EE2 could certainly result in 
population collapse due to a complete loss of male brood- 
ing in pipefish. However, under the lowest levels of EE2 
exposure investigated in the present study, males were able 
to become pregnant and the nature of sexual selection act- 
ing on females remained mostly unchanged relative to con- 
trol populations. 

The most important observation in our study indicated 
that EE2 exposure affected the reproductive success of 
females, but the proximate effect appeared to be positive 
rather than the negative effect we would have predicted a 
priori. Females in the control group that mated multiply 
transferred fewer eggs per mate than females with only one 
mate, indicating that females may have rationed their eggs 
or become egg limited as they mated with additional males. 
In the case of EE2-exposed females, however, we found a 
different pattern: multiply mated females transferred a 
comparable number of eggs per mate compared with singly 
mated females. When we look at the total number of eggs 
produced by females, we see that females exposed to EE2 
that mated multiply produced many more eggs overall than 
the multiply mated females in the control group. The mul- 
tiply mated females did not differ in body length between 
the treatments, so this difference cannot be attributed sim- 
ply to larger females having higher fecundity, a trend typi- 
cally seen in other species of pipefish (Braga Goncalves 
et al. 2011). The most likely explanation for this observa- 
tion is that exposure to low levels of EE2 stimulates egg 
production in female Gulf pipefish thereby increasing their 
reproductive rates and that females can only realize these 
enhanced reproductive rates by mating with multiple 
males. This interpretation is consistent with observations 
involving fathead minnows, in which short periods of EE2 
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exposure have been shown to increase egg production 
(Jobling et al. 2003). However, chronic low levels of expo- 
sure eventually cause decreased egg production in fathead 
minnows (Jobling et al. 2003). Several studies have also 
shown reproductive failure in zebrafish as a result of either 
long-term exposure to low levels of EE2 or brief exposure 
to higher levels (Van den Belt et al. 2003; Nash et al. 2004; 
Xu et al. 2008), raising the possibility that the short-term 
benefit that we observed in pipefish could be offset by 
reduced lifetime fitness for EE2-exposed females. Future 
studies involving a longer time frame will be necessary to 
address this question. 

This increase in reproductive rates in EE2-exposed 
females also appeared to have an effect on the trade-off 
between pre- and postmating episodes of selection. In the 
control group, females that mated multiply were able to 
transfer more eggs than singly mated females over the 
entire experiment but in turn had more of their eggs fail to 
develop. This trade-off between greater mating success in 
the premating phase of selection and a decrease in offspring 
survivorship, a postmating mechanism, was not present in 
the EE2 multiply mated females. This observation raises the 
possibility that EE2-exposed females produce eggs of higher 
quality as a result of the excess estrogen resource. It is also 
important to note that while we did not dissect out the 
ovaries of unmated females, there did not visually appear 
to be any degradation of functioning ovaries in the females 
that did not mate. 

As noted above, while an exposure of 2 ng/L seems to 
benefit females in the short term, EE2 exposure at 5 ng/L, a 
concentration on the lower side of the range of EE2 
detected in the natural environment, results in devastating 
reproductive impacts on Gulf pipefish populations. 
Females are largely unaffected, but males are seriously com- 
promised by this level of EE2 contamination. In our pilot 
study, males exposed to 5 ng/L failed to carry pregnancies 
to term, were unable to mate, and showed abnormal brood 
pouch morphology. At higher levels of exposure, the effects 
are even more dramatic. For instance, male Gulf pipefish 
exposed to 100 ng/L displayed secondary sexual traits that 
are normally only found in females, including iridescent 
bands and a deeply keeled abdomen. Furthermore, these 
males showed a reduced ability to mate even after being 
removed from the short-term EE2 exposure, with a mini- 
mum lag time of 4 days until pregnancy (Partridge et al. 
2010). Thus, we predict that exposure in natural popula- 
tions approaching 5 ng/L or higher will reduce the repro- 
ductive potential of Gulf pipefish populations and that the 
negative effects will be mediated almost entirely by the 
impact of EE2 exposure on males rather than on females. 
Several studies have documented levels of EE2 ranging 
from 5 ng/L to below, indicating that both of the concen- 
trations in our study are within the range of possible EE2 
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contamination occurring in the environment where natural 
populations of the Gulf pipefish reside. 

We did not see a breakdown in selection in the Gulf 
pipefish's sex-role-reversed mating system in populations 
exposed to EE2. Instead, we see no significant changes in 
the opportunity for selection in both male and female Gulf 
pipefish. As a result of increased female reproductive rates 
in EE2, we actually see a small but nonsignificant increase 
in selection acting on body length in both males and 
females. In a previous study by Saaristo et al. (2009b) doc- 
umenting a breakdown in sexual selection in sand gobies, a 
species with strong sexual selection on males, sexual selec- 
tion was disrupted due to the feminization of males. If the 
male traits that sand goby females use to choose their 
mates, such as aggressiveness and courtship displays, are 
disrupted by EE2 exposure, then sexual selection would be 
expected to break down. In our sex-role-reversed system, 
we did not detect a breakdown of selection as a result of 
male feminization. If anything, selection on females was 
slightly stronger in the EE2 treatment than in the control 
(as evidenced by the significantly steeper Bateman gradient 
for EE2-exposed females), which could be due to the 
greater number of eggs available and their increased survi- 
vorship, possibly causing an increase in female-female 
competition. Sand gobies have conventional sex roles with 
male-male competition for access to females, so it makes 
sense that feminization of males would reduce their ability 
to compete for mates. Part of the reason we found a differ- 
ent effect of EE2, compared with the complete collapse of 
selection acting on male sand gobies, certainly stems from 
the sex-role-reversed mating system of Gulf pipefish. 
Increased feminization of females appears to have made 
them more fecund and more able to compete for mates. 
However, feminization of male Gulf pipefish did not seem 
to affect the mating system until the males began to lose 
their ability to maintain functional brood pouches, a situa- 
tion that results in a complete cessation of reproduction 
rather than quantitative changes in the intensity of sexual 
selection. It is also important to note that in the study con- 
ducted by Saaristo et al. (2009b), female sand gobies were 
not exposed to EE2. However, in the present study, we 
exposed both sexes to EE2 to investigate the effects of expo- 
sure in a setting that would more closely mimic contamina- 
tion of the natural environment. 

In summary, low levels of EE2 exposure enhanced repro- 
duction in female Gulf pipefish and increased the fitness of 
multiply mated females during both pre- and postmating 
episodes of selection. EE2 exposure did not disrupt premat- 
ing sexual selection in the Gulf pipefish mating system; if 
anything, sexual selection was slightly stronger in the 
exposed populations, owing to the increased fecundity of 
exposed females. Even though this study documented the 
effects of low levels of EE2 exposure on the Gulf pipefish 
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mating system, it is critical to put the level of exposure in 
perspective. In our study, EE2 had a positive effect on 
female egg production at low levels of 2 ng/L. However, at 
EE2 exposure levels of 5 ng/L and higher, male receptivity 
decreased and male fitness plummeted to zero. In conclu- 
sion, this study has added to our knowledge on the effects 
of EE2 on pipefish, providing more insight into the effects 
EE2 has at the population level rather than on the individ- 
ual level. Now that we have established that successful mat- 
ings can occur under low levels of EE2 exposure, the next 
question is how these pollutants affect these fish at various 
stages of their life cycle, especially in terms of juvenile 
recruitment and population viability in EE2-contaminated 
waters. 
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